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ABSTRACT 



A five- level baroclinic model for numerical prediction of height 
fields at 1000, 850, 700, 500, and 300 mb is developed and programmed. 
Several forms of the model are tested to determine the contribution 
of the various terms of the vorticity equation. Forecasts from three 
different methods are compared. The model shows utility in prognosis. 
Although there are individual differences among the methods, no one of 
the three is significantly better over the entire hemisphere on the 
data tested here. 

The writers wish to express their appreciation to Professor George 
Haltiner of the U. S. Naval Postgraduate School for his assistance and 
guidance in this investigation. Appreciation is also expressed to the 
personnel of the U. S. Navy Fleet Numerical Weather Facility for their 
cooperation during the preparation of this paper. We are especially 
indebted to CDR L. C. Clarke and LT G. Lawniczak for use of their 
solution of the omega equation, and for their assistance to us, and 
to LCDR M. Jane Frawley for her aid and advice in programming. 
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Table of Symbols 



the upward component of the apparent gravitational acceleration 

the altimeter correction P -'Z-'Z-y 

the height of an isobaric surface 

the pressure altitude in the standard atmosphere 

the Coriolis parameter, 3 _Q where <j is the geograph- 

ical latitude 

the mean value of the Coriolis parameter, fl Sm H 5 

the geopotential 

the relative vorticity £ = ^ 

the absolute vorticity £ + 'f 

the stream function for the non-divergent component of velocity 

the velocity potential for the divergent component of velocity 

the vertical velocity of air in pressure coordinates 

the atmospheric pressure in millibars 

T 50 

the static stability <T - “ 
the universal gas constant 
grid distance, 381 km at 60N 
map scale factor 

del operator on a constant-pressure surface 
horizontal Laplacian operator on a constant-pressure surface 
horizontal Jacobian operator T(A> b') 2 ^ ^ y— 

finite difference Laplacian 



V 2 A=(A l4 ,.+A (j j +l tAi.g+Ag, l -4A, :j ) 



finite difference Jacobian 






1. Introduction. 



Until very recently, the barotropic models have been the principal 
source of numerical weather prognoses. Their shortcomings in areas of 
baroclinic development, plus the advent of sophisticated digital computer 
systems, have spurred interest in baroclinic models. It is the purpose 
of this research to develop and program a five-level baroclinic model and 
to investigate the contribution of the various terms of the vorticity 
equation to the prognosis of the pressure field. 



1 



2. Background* 

In the model developed here, prognosis is accomplished by means of 
the vorticity equation, which, with pressure as the vertical coordinate, 
can be written 
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Basic to the derivation of this equation is the Helmholtz theorem which 
expresses the wind as the sum of a divergent and a non- divergent compon- 
ent: 

V = IK * \7 x + \/ x. . ( 2 ) 

The potential function X for the divergent part of the wind is obtained 
through solution of the continuity equation 

V a X + | *« o. a, 

The stream function ~ K j r for the non-divergent portion is approximated by 
solution of the balance equation 



f v 'y+4fX!X(-!X)l+ vf-vf- 



(4) 



as described by G. 'Arnason [1]> where (Q is the geopotential. Use of 
the stream function in the vorticity equation enables one to avoid the 
geostrophic approximation and the concomitant spurious anticyclogenesis . 

The vertical component of relative vorticity, £ , is the Laplacian 
of the stream function 

c-xr (5 
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and absolute vortlcity, '/ , Is the sum of this and the vertical component 
of the earth's vorticity, f. 

v = r + f (6) 

Vertical velocity, 60, is computed by the diagnostic equation 



wt r« a P=- f 






(7) 



developed by Clarke and Lawniczak [3] . The forcing function is calculated 
with the geostrophic approximation. This results in some inconsistencies 
in the model but to avoid reprogramming the entire equation, we have gone 
ahead with comparison- type experiments. Vertical velocity is assumed zero 
at 100 mb, but other boundary conditions are among the options available 
to the user of this computer program as follows: 

a) the lower boundary can be level or include terrain; 

b) vertical velocity can be set equal to zero at the lower boundary, 
or computed with or without terrain effects and with or without the 
effect of friction; 

c) the stability parameter ( T can be a point variable, a function of 
pressure alone, or a constant; and 

d) the maximum acceptable residual in the 3-dimensional relaxation 
can be selected by the user. 

In the vorticity equation, the contributions from the twisting term, 

f \ 

[ TIT c>x 2 p — TV d y dp ) an< * t ' ie vertical advection of vorticity term, 

00 ^ are usualL y considered small, and some investigations have 

shown that they tend to compensate each other [2]. Also, integral theo- 
rems show [4] that to avoid the mean generation of vorticity in a closed 
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volume, the two terms should be included or neglected together. Both 
have been omitted from the initial investigation reported here. 

A further simplification involves omission of the advection of abso- 
lute vorticity by the divergent portion of the wind, yy'- + " r y" • 

For consistency, it is necessary then also to substitute the mean value 
of the Coriolis parameter, f , for the absolute vorticity in the diver- 
gence term, f . 

In the equations that follow, the stream function has been redefined 
so as to be in units of height, 

y* =* j-E . (8> 

The resulting forms of the vorticity equation, expressed in finite- 
difference notation, are 



/wo* 9 vt - 7 1 ^^ 
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and the complete equation 



rrrf 0 rpt AIL- _ - Cj/yrf TT /-. I I 9^ 

ol 1 ^ Y At f 4 cl ‘ J/VP.V; [ Ux Ax Ay AY/ ‘ A p 

9 /W^ ptp^E. , g /W>Y A^ \ 

f d 1 AP f i Uy AYAp A* AyAp/ ; ' 

where d is the grid distance and /W\ is the ratio of this distance to the 
corresponding geographical distance. 
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Horizontal spatial operators are approximated by the usual centered 



differences • To approximate the vertical derivatives in the vorticity 
equation a Lagrangian parabola was fitted through three points. For a 
centered difference the equation is 



All computations were done on the Control Data Corporation model 
1604 computer, using the U. S. Navy Fleet Numerical Weather Facility 
octagonal grid. This grid, centered at the pole and inscribed in the ION 
latitude circle, consists of 1977 grid points with a grid size of 381 km 
at 60N. 

The week of 22 November 1962 was a period of considerable cyclogene- 
sis in dense-data areas, and therefore excellent for testing a baroclinic 
model. Initial data used were the U. S. Navy Fleet Numerical Weather 
Facility height analyses at 1000, 850, 700, 500, and 300 mb for 1200Z on 
22 November 1962. The height analyses for 0000Z and 1200Z 23 November 
were used in the evaluation of the predictions. 
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3. 



Procedures . 



In the diagnostic 60 -equation for the initial tests discussed here, 
the boundary conditions used are: 

a) vertical velocity equal to zero at 100 mb; 

b) vertical velocity equal to zero at the lower boundary; 

c) no terrain effects at the lower boundary; and 

d) the stability parameter, <j~ , is a function of pressure. The pro- 
gram is used in n forecast mode"; that is, temperatures are computed hydro- 
statically, rather than taken from objective analysis. The maximum 
residual in the 3-dimensional relaxation is 1 x 10”^ mb per second. 

The sequence of calculations programmed can best be explained with 
reference to figure 1. Input data in every case are the U. S. Navy Fleet 
Numerical Weather Facility height analyses at the five standard levels, 
modified to ensure vertical consistency. 

In the first method with equation (9), stream functions are used 
throughout. The inversion to height fields [l] , effectively the reverse 
of the balance equation, is thus bypassed and stream functions are used 
for input data for both diagnostic and prognostic equation. This is 
hereafter referred to as method I. The resulting 24-hour prognoses are 
displayed on charts 6 through 10 and discussed in section 4 below. 

Since the geostrophic approximation is used in the forcing function 
of the ^-equation, and temperatures are computed hydrostatically, there 
might be some reason to suspect the vertical velocities computed from 
height fields could be more accurate than those from stream functions, A 
completely geostrophic model was programmed, using pressure-height data 
throughout. This method showed the expected spurious anticyclogenesis 
even at 12 hours and was carried no further. 
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The complete process diagrammed in figure 1 has been used in the re- 



mainder of this investigation. In the first time step the initial data 
(five vertically-consistent height fields) are used as input to the bal- 
ance equation and to the CO -equation. From the former are obtained stream 
functions for the five levels; the latter yields vertical velocities at 
the same levels. With these stream functions and vertical velocities the 

vorticity equation (equation (9) for method II, equation (10) for method 
E 

III) gives , the local change in stream function, at the five levels* 

At 

When these are added to the initial values, the first forecast time step 
is complete. As shown in figure 1, these new E values are now input 
fields for both the vorticity equation and the balance equation; the 
latter being used for inversion to pressure height as required by the 
CO -equation. 

Subsequent time steps follow the same routine, except in the applica- 
tion of A E ^ which is done according to the centered-dif ference equation 



In order to minimize the errors inherent in the initial forward 
time step, the first time increment is only 15 minutes. Centered time 
steps are then used to complete successive 15-minute prognoses to one 
hour, and finally 1-hour prognoses to 12 hours. The fields are then 
smoothed lightly to remove any small perturbations resulting from trunca- 
tion or round-off errors, and the process is repeated to 24 hours* 



At 




(13) 
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4. 



Results . 



In review: method I uses stream functions throughout, and with 

equation (9) as the prognostic equation excludes the advection of abso- 
lute vorticity by the divergent portion of the wind and substitutes f 



fields to pressure heights at every time step, but retains the other 
simplifications of method I; method III also inverts to pressure height 
at each time step, but uses the more complete equation (10) as the prog- 
nostic equation. 

Charts 1-5 are the pressure height analyses at the beginning of the 
period. The next 15 charts are the results of the three 24-hour fore- 
casts: method I in charts 6-10; method II, charts 11-15; and method III, 

charts 16-20. Charts 21-25 are the observed height analyses at the end 
of the forecast period. Chart 26 is the U. S. Navy Fleet Numerical 
Weather Facility barotropic forecast at 500 mb for the same period. 

Table 1 shows the contour origin and interval for all the charts. 

In general, corresponding charts from the three methods are quite 
similar. They all retain a typical meteorological appearance and the 
various systems show a reasonable vertical structure. Further indica- 
tion of their similarity is contained in table 2. The ’’pillow” referred 
to is 



and the RMSE is the root-mean- square error calculated in the usual 
fashion: 



for ?l in the divergence term; method II inverts the stream- function 



Pillow 




( 16 ) 



n 7 7 



8 



RMSE 
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Nevertheless, a level-by-level comparison of the three forecast 
methods shows some differences among them. At 1000 mb, for instance, the 
two low centers over North America merged and deepened, while the trough 
moved nearly 20° east. All methods forecast the merging and consider- 
able trough movement, with method I giving the best placement of the 
trough. During the forecast period a low center developed in the area 
of weak gradient southwest of the Kamchatka peninsula. This development 
was forecast by all methods, although with too much cyclogenesis. 

Method I was again the best. In no case is there any evidence of spuri- 
ous anticyclogenesis, even at low latitudes. 

The forecasts at 850 and 700 mb show movement consistent with that 
forecast at 1000 mb. By 700 mb the tilt of the baroclinic systems has 
become apparent. A good example is the stationary low center over the 
Northwest Territories, which deepened by 26 meters. All forecasts held 
it stationary, but methods I and III forecast two and seven meters deep- 
ening, respectively, while method II forecast a 2-meter filling. The 
weakening of the low center east of Greenland was forecast by all methods, 
but best by method II. 

At 500 mb, since the chart of the U. S. Navy Fleet Numerical Weather 
Facility barotropic forecast has been included, there are four charts to 
compare. There is a striking similarity between the barotropic progno- 
sis and that from method I. The method II and III prognoses, themselves 
much alike, are quite different from the other two. Method I and the 
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USNFNWF model both forecast the movement and reorientation of the trough 
over Eastern United States, and the weakening and displacement of the 
low center near Hudson Bay. Methods II and III show the trough lagging 
far behind its true position, and retain the low as a stationary system.. 
Method III even deepened it slightly. In contrast, methods II and III 
produced much better forecasts than the others in the western Europe- 
Mediterranean area. 

At 300 mb all three methods agree in general on the movement and 
development of the major systems. This development is in some cases out 
of phase with the actual development, but no one method can be selected 
as better than the rest. The initial and verifying charts at 300 mb 
show more effect from the process of making the data vertically consist- 
ent than was evident at lower levels. 

Referring again to table 2, there seems to be no best method among 
the three tested. Certainly the addition of the advection of absolute 
vorticity by the divergent wind and the use of absolute vorticity in the 
divergence term produced little change and no over-all improvement in 
this particular forecast. The effect of the inclusion of the vertical 
advection of absolute vorticity and the twisting term remains to be 
tested. 
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Figure 



Table 1 



Level 


Contour origin 


Contour interval 


(mb) 


(meters) 


(meters) 


1000 


120 


60 


850 


1440 


60 


700 


3060 


60 


500 


5580 


60 


300 


9120 


120 
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Table 2 



Method 


Level 

(mb) 


Pillow 

(meters) 


RMSE 

(meters) 


i 


1000 


- 7.9 


43.6 




850 


- 1.8 


37.6 




700 


- 1.8 


41.8 




500 


- 1.2 


51.3 




300 


13.4 


80.3 


ii 


1000 


-11.3 


43.3 




850 


- 2.7 


40.8 




700 


0 


42.7 




500 


- .6 


53.4 




300 


11.9 


77.3 


in 


1000 


- 9.1 


49.1 




850 


- 3.4 


43.0 




700 


- 2.4 


43.6 




500 


- 1.5 


54.4 




300 


13.4 


80.0 


USNFNWF 


500 


- 4.9 


47.0 


Persistence 


500 


- 3.3 


57.4 
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APPENDIX I 



The following charts were computed during the course of this study. 
They are included here for easy reference. 



Charts 1-5 


Analyses, 12Z 22 Nov 1962 


Charts 6-10 


24-hour prognoses, method I 


Charts 11-15 


24-hour prognoses, method II 


Charts 16-20 


24-hour prognoses, method III 


Charts 21-25 


Analyses, 12Z 23 Nov 1962 


Chart 26 


24-hour prognosis, USNFNWF 
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